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Abstract 

O 

By observing reconstituted chromatin by fluorescence microscopy (FM) and atomic force mi- 

^? 

croscopy (AFM), we found that the density of nucleosomes exhibits a bimodal profile, i.e., there is 

> 



a large transition between the dense and dispersed states in reconstituted chromatin. Based on an 
analysis of the spatial distribution of nucleosome cores, we deduced an effective thermodynamic 
potential as a function of the nucleosome-nucleosome distance. This enabled us to interpret the 
folding transition of chromatin in terms of a first-order phase transition. This mechanism for the 
condensation of chromatin is discussed in terms of its biological significance. 

PACS numbers: 87.15.He, 64.75.+g, 36.20.Ey 
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Genomic DNA in eukaryotes is compactly folded into chromatin through several hierar- 



chical packings 



lj. The fundamental unit of such packing, the nucleosome, consists of 146 bp 



of DNA wrapped around a histone octamer (two molecules each of H2A, H2B, H3 and H4). 
The semi- flexible DNA chain wraps around a histone core about two turns |2j. It is widely 
expected that the manner of packing and the dynamics of nucleosomes are associated with 
gene activities in living cells [3-6]. There have been many studies on the static [7-10] and 



dynamic 



ic ll| 



properties of nucleosomes. It has been shown that nucleosomes condense under 



121. However, the nature of 



various conditions; for example, under a high salt concentration 
the higher-order structure is not well understood and the underlying physics of nucleosome 
condensation have not yet been clarified. It has been suggested, for example, that mea- 
surement of the actual interaction energy between nucleosomes is essential [8,11,13,14] for 
obtaining deeper insight into chromatin condensation, but there have been no experimental 
studies on the interaction potential. 

On the other hand, it has recently been found that linear DNA larger than several tens 
of base pairs exhibits a large discrete coil-globule transition, accompanied by a change in 
density on the order of 10 4 -10 5 , upon the addition of various kinds of condensing agents 
[15,16]. It has been revealed that this transition is a first-order phase transition under the 
criterion of Landau, i.e., an ON/OFF transition from an elongated coil state to a compact 



globule state 



171 ] . Very recently, such an ON/OFF transition of DNA has been suggested 



the coil-globule transition is 



id- 



to play an important role in gene activity 17]. Discreteness o: 
a general characteristic of single semi-flexible polymer chains 

The purpose of this study was to obtain a deeper understanding of the conformational 
changes of chromatin, which is essential for obtaining insight into genetic activity including 



duplication, transcription, etc. We investigated the physical properties of reconstituted chro- 



matin using fluorescence microscopy (FM) and atomic force microscopy (AFM). We obtained 
the pair interaction potential between nucleosome cores from an analysis of AFM images 
and used it to dissect the mechanism of chromatin compaction, or chromatin condensation 



191. 



The preparation of core histones and 106-kbp plasmids (circular DNA) and the reconsti 



tution of chromatin were carried out as previously reported 



In this study, the mass ratio 



[histone]/[DNA] was set to 1.0 and 1.3. Concentration of NaCl, which affects the interac- 
tion between nucleosomes 7[, was 50 mM throughout this article. Reconstituted chromatin 
samples were fixed with 0.3% glutaraldehyde in lOmM Hepes-NaOH [pH 7.5] for 30 minutes 
at 25 ° C. The chromatin was placed on a freshly cleaved thin piece of mica (thickness; ca. 
30-50 fim) stuck to a glass cover plate (Matsunami Glass, No. 1, Japan) for FM obser- 
vation. The mica surface was pre-treated with 10 mM spermidine. This treatment was 
performed so that chromatin would adsorb onto the mica surface. The chromatin complexes 
were visualized by fluorescent microscopy using a fluorescent dye, 0.1 /iM 4', 6-diamidino- 
2-phenylindole (DAPI). The observation was performed on a droplet (10 /xL) instilled on 
mica. The sample droplet on mica was washed with Millipore water and blown dry with 
nitrogen gas for 5 minutes. Fluorescent chromatin complexes were observed under a Zeiss 
Axiovert 200 microscope with a 100 x oil-immersed objective lens at 25 ° C, and recorded 
on Axio Vision with an AxioCam camera. To obtain two-dimensional real-time fluorescent 
image data, an inverted microscope (IX-70, Olympus) with a lOOx oil-immersed objective 
lens and a highly sensitive EB-CCD camera with an image-processing system (Hamamatsu 
Photonics) were used. The video data were recorded on videotapes, and then analyzed with 
personal computers. Due to the blurring effect [21[ in the observation with a highly sensitive 
video system, the size of observed DNA images was assumed to be slightly larger (ca. 0.3 



/iin) than the actual size of the chromatin. The reconstituted chromatin structures observed 
by FM observation were analyzed under Tapping Mode™ in AFM (Nanoscope Bioscope, 
Digital Instruments) in air at room temperature. Both FM and AFM images were obtained 
on exactly the same chromatin. The distance distribution of nucleosomes was obtained on 
AFM images acquired with a high- resolution AFM apparatus (Nanoscope Ilia). 

Figure 1 shows quasi-three-dimensional (3D) FM images (A-C), corresponding 2D images 
(insets), and AFM images (D-F) of reconstituted chromatin adsorbed onto a mica surface. 



221. The 



The pictures in Fig. 1 (A, D), (B, E), and (C, F) show exactly the same molecules 
mass ratio [histone]/[DNA] is 1.0 in (A, D) and 1.3 in (B, C, E, and F). In (D), nucleosomes 
are dispersed in the chromatin, whereas in (E) the condensed and dispersed parts coexist. 
In (F) the chromatin is entirely condensed. Under the condition [histone]/[DNA] = 1.3, the 
partially and entirely condensed states exist in almost equal proportions, while there is a 
very low proportion (< 5%) of the dispersed state. Under [histone]/[DNA] = 1.0, more than 
95% of chromatin is in the dispersed state, while the remainder is in the partially condensed 
state. Despite the low resolution, FM observation provides information on the degree of 
condensation of individual reconstituted chromatin in the bulk solution as well as on the 
surface. 

To evaluate the actual size of reconstituted chromatin in bulk solution, we measured the 
Brownian motion of individual chromatin complexes using FM. From the time-dependence 
of the mean square displacement of the center of mass of chromatin, we obtained the dif- 



fusion constant D using the following relationship: ((r(0) — r(t)) 2 ) = 4Dt 23j. The hy- 
drodynamic radius Rh is thus deduced as in Table I using the Stokes-Einstein relationship, 
Rh = kBT/(6irr]D), where ks is the Boltzmann constant, T is the absolute temperature 
and 7] is the viscosity of the solvent (0.89 mPas for pure water at T = 297 K). We also 



evaluated the size of reconstituted chromatin by AFM measurement. The major axis Rl 
and the minor axis Rs of reconstituted chromatin were measured in two-dimensional AFM 
images. Table I shows Rl and Rs, together with -Rafm = VRlRs- The relative ratio of 
Rafm/Rh is almost the same for [histone]/[DNA] ratios both of 1.0 and 1.3, regardless of 
the difference in the observed state, i.e., the former is chromatin on a surface and the latter 
is chromatin in bulk solution. 

We obtained the distribution of distances N(r) as a function of the distance r between 
pairs of nucleosomes in AFM images. It is expected that nucleosomes exhibit a nearly 



l — l 



equilibrium structure on a 2D mica surface 



24j . We assumed that the radial distribu- 



tion function g(r) is proportional to N(r)/r. Since the nucleosome density is low enough 
with [histone]/[DNA] = 1.0 (see Fig. 1(A)), the pair potential can be roughly deduced 
from g(r) assuming a Boltzmann distribution, U(r) = — /c^Tlng^r) [25-27], where g(r) and 
U(r) were calculated after fitting N(r) to a polynomial. The analysis was performed for 
[histone]/[DNA] = 1.0, since in this condition we can count all of the nucleosomes in recon- 
stituted chromatin. The result is shown in Fig. 2 and indicates that U (r) has a minimum 
0.3 ksT depth at r pa 13 nm. By fitting the experimental data in Fig. 2 with the following 
equation of potential, U(r) = —Ar~^ + Br~ u (/x, v: integer; /i < v), we obtain A = 2.4, B 
= 12, // = 9, and v = 23, where r is normalized with the diameter of a nucleosome (11 nm). 
These values are much larger than those of the Lennard- Jones potential (// = 6 and v = 
12). The large values of the exponents can be attributed to the large excluded volume of the 
nucleosome core. A similar profile of potential, i.e., relatively large values for the exponents 
\i and v, can also be deduced from an analysis of the interaction of neighboring nucleosomes 
along the DNA chain (data not shown). Such large values of the exponents can explain the 
discrete nature of the transition as described below. 
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Next, we will discuss the conformational stability of chromatin based on the above pair 
interaction energy. The total free energy F of reconstituted chromatin with n nucleosomes 
can be described as 

F{n) = F dA + F ist (n), (1) 

where F e \ a is the entropic elasticity of the DNA chain and Fi nt (n) is the volume interaction 
between nucleosomes. We neglect the volume interaction of double-stranded DNA, since the 
thickness of DNA (~ 2 nm) is much smaller than the diameter of a nucleosome (~ 11 nm). 
Using the swelling parameter a, we obtain 

F ei Jk B T = a 2 + a~ 2 , (2) 

where a 2 = (R 2 )/(Rl). (R 2 ) is the mean square of the radius of gyration and R is the 
analogous size of an ideal Gaussian coil, a 2 and a~ 2 correspond to extension and compres- 
sion of the chain, respectively |28(. Assuming that U(r) has a narrow minimum so that 
nucleosome-nucleosome interaction occurs only among the nearest neighbors, we can take 
Fiut(n) ~ nU(r). By adapting U(r) = —Ar~ 9 + Br~ 23 , Fi nt (n) can be written as 

F iDt (n, p)/k B T ~ -Anp 3 + Bnp 7 - 7 , (3) 

where p is the normalized density of nucleosomes, p ~ r~ 3 . Thus, we obtain 

F(n, p) /k B T ~ a 2 + a" 2 - Anp 3 + Bnp 7 - 7 . (4) 

In the present model, Rq decreases with an increase in n. Thus, we define L as the apparent 
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contour length of chromatin, L = L — an, where L is the contour length of DNA without 
any histones (36 pm for this sample), and a is the length of DNA wrapped around a histone 
octamer, ca. 50 nm (146 bp). With the Kuhn length A (100 nm for a DNA chain) and the 
number of Kuhn segments N s , the size of an ideal Gaussian chain is described as (i? 2 ,) 1 / 2 = 
XNl /2 . Since N s = L/X, we obtain a modified R , which we call R' , {R'£} 1/2 = A(L/A) 1/2 = 
\ l l 2 L 1 / 2 . Using the relation p ~ ni? -3 , we obtain 

HTT = V3»r , + + (5) 

nk B T n L ' d X{L — an) n b ' 6 

It is obvious that F(n, pj/nksT has two minima. The one at the lower-density region is 
derived from the first and second terms in Eq. (5) and the other at the higher- density region 
is derived from the third and fourth terms. As n increases, the minimum of the condensed 
state becomes deeper while a double-minimum profile is maintained. When the two minima 
have a similar depth, two different states of high and low nucleosome density coexist. Figure 
3 shows the free-energy profiles of nucleosomes with n = 400, 500 and 600 calculated with 
Eq. (5) together with schemes of the corresponding conformations in three dimensions. 
The observed result regarding the elongated conformation in Fig. 1 corresponds to the 
deeper minimum at p ~ 0. With an increase in n, the free energy of the condensed state 
becomes the absolute minimum, which reproduces the experimental trend. In the actual 
experiment for [histone] /[DNA] = 1.3, there are two different states: fully condensed and 
partially condensed. Any instability due to the interfacial energy between the condensed and 
dispersed parts should be negligible considering the zero-dimensional nature of the interface. 
Therefore, we can expect the appearance of an intrachain phase-segregated state as in Fig. 
1(B). Previous studies on native linear DNA chains have indicated that a phase-segregated 
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state is actually observed in individual DNA molecules and that the characteristic scale of 
segregation depends on the degree of the surviving electronic charge in the condensed part 
[29,30]. 

Several biological studies have demonstrated that the higher-order folding of chromatin 
fiber and its dynamic structural changes largely depend on the proper functions of various 
structural and regulatory proteins in the nucleus [31,32], which are critical for gene expression 
and chromosome segregation. On the other hand, the results obtained in this study, together 
with those in several previous studies [11,33], suggest that the physical properties of a DNA 
strand (length and superhelicity) and the interaction between nucleosomes play fundamental 
roles in chromatin dynamics. The higher-order architecture of chromatin is determined by 
the fundamental properties of chromatin fiber itself. In this sense, it should be noted that 
chromosomes are composed of several chromatin loops on the order of ~100 kb [34-36], 
which is approximately the same length as used in this study. A first-order large-scale 
conformational transition may explain why previous experiments have failed to observe 
an intermediate state in chromatin condensation. More importantly, it may explain the 
switching of a large number of genes. Recently, it has been reported that transcription 
is completely inhibited through an all-or-none transition in the structure of giant DNA 
molecules [17,37]. Thus, we would like to propose a hypothesis of gene self-regulation through 
a combination of large-scale ON/OFF switching due to a conformational transition and 
several regulatory factors to bind specific regions. 
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versity and Universality in Physics", the Special Co-ordination Funds (to K. T.), and the 
COE Research Grant from Ministry of Education, Culture, Sports, Science, and Technology 
of Japan (to K. T.). S. H. Y was the recipient of a predoctoral fellowship from the Japan 
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TABLE I: The major axis of a chromatin R^, the minor axis Rs, Rafm, an d the hydrodynamic 
radius Rh- Rl and Rs were measured in AFM images by assuming that a chromatin complex 



is elliptical. Rh was measured in fluorescent images. -Rafm ( = V^lRs ) an d Rh are the 
hypothetical radii of chromatin if they are considered to be spheres on a mica surface and in bulk 
solution, respectively. 



[histone]/[DNA 


Rl (nm) 


R s (nm) 


-Rafm (nm) 


Rh (nm) 


1.0 


420 ± 180 


250 ± 90 


320 ± 120 


370 ± 70 


1.3 


220 ± 100 


130 ± 90 


170 ± 100 


190 ± 50 
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Figure Captions 

FIG. 1: Reconstituted chromatin with [histone]/[DNA] = 1 (A, D) and 1.3 (B, C, E, and F). 
The images in (A, D), (B, E), and (C, F) are photographs of the same chromatin complex. 
(A-C) Light-intensity distribution of fluorescent images of reconstituted chromatin situated 
on a mica surface. Insets are corresponding 2D fluorescent images. (D-F) AFM images of 
the same chromatin as in the fluorescent images. The scale bar is 0.2 /iin in AFM images. 

FIG. 2: The pair potential of nucleosomes obtained by an equation for the Boltzmann 
distribution. Inset: The nucleosome-nucleosome distance distribution N(r) of chromatin 
reconstituted with [histone]/[DNA] = 1.0. U(r) is deduced by polynomial curve-fitting of 
N(r). 

FIG. 3: Free-energy profiles of a nucleosome with n = 400, 500, and 600 and corresponding 
schematic representations, where n is the number of nucleosomes in a single chromatin 
complex. The free energy as a function of the normalized density of nucleosomes in the 
three-dimensional conformation was calculated from Eq. (5). 
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